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An effect of nondissipative drag of a superfiuid flow in a system of two Bose gases confined in two 
parallel quasi two-dimensional traps is studied. Using an approach based on introduction of density 
and phase operators we compute the drag current at zero and finite temperatures for arbitrary ratio 
of densities of the particles in the adjacent layers. We demonstrate that in a system of two ring-shape 
traps the "drag force" influences on the drag trap in the same way as an external magnetic flux 
influences on a superconducting ring. It allows to use the drag effect to control persistent current 
states in superfluids and opens a possibility for implementing a Bose analog of the superconducting 
Josephson flux qubit. 
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I. INTRODUCTION 

The existence of nondissipative supercurrents is a common feature of superconducting and superfluid systems. 
Among various applications of superconductivity considerable attention was given to the use of superconducting 
circuits as very sensitive magnetometers (superconducting quantum interferometer devices). At present an interest 
to such systems is renewed in view of a possibility to use superconducting circuits with weak links as elements of 
quantum computers (Josephson qubits). Basing on a similarity between superfluids and superconductors one can 
expect that the former ones may also be used for implementing qubits. 

Supercurrent in superconductors is coupled to vector potential of electromagnetic fields. It allows to control per- 
sistent current states by external fields. Obviously, there is no such a channel for a control in uncharged superfluid 
Bose systems. In this paper we study another possibility based on a nondissipative drag effect. 

The drag in normal systems has been investigated experimentally and theoretically by many authors (see, for 
instance, reviews 0, 0). The main attention was given to the study of bilayer electron systems in semiconductor 
heterostructures. In such systems an interlayer drag effect takes place. The effect is caused by electron-electron 
scattering processes and it reveals itself in an appearance of a drag voltage in one layer when a normal current flows 
in the adjacent layer. If the former layer is in a closed circuit, the drag voltage induces a drag current flowing through 
the circuit. The effect is accompanied by a dissipation of energy and takes place only at finite temperatures. Roughly, 
the drag voltage increases by the T 2 -law (deviation from this law, observed experimentally is connected with a 
phonon contribution to the interaction between the carriers 0). 

In superfluid and superconducting systems another kind of drag may take place. This drag is nondissipative and is 
connected with a redistribution of a supercurrent between two superfluid (superconducting) components. In difference 
with the drag in a normal state the superfluid drag has the largest value at zero temperatures and decreases under 
increasing the temperature. On the existence of nondissipative drag in superfluid systems was pointed out for the 
first time in the paper by Andreev and Bashkin [f| . In this paper a three- velocity hydrodynamic model of a 3 He- 4 He 
mixture was developed. It has been shown that superfluid behavior of such systems can be described under accounting 
the "drag" term in the free energy. This term is proportional to the scalar product of superfluid velocities of two 



components times the difference between the effective and the bare masses of 3 He atoms. In the paper by Duan and 
Yip the nondissipative drag effect in superconductors has been studied. The authors of [|| argue that the value 
of drag can be obtained from the energy of zero-point fluctuations. It was shown that this energy contains a "drag" 
term analogous to one obtained in [5| in the hydrodynamic approach. The theory of the nondissipative drag in a 
bilayer system of charged bosons was developed by Tanatar and Das Q and by Terentjev and Shevchenko §. 
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The existence of nondissipative drag in a system of two one-dimensional wires in a persistent current state was 
predicted by Rojo and Mahan . It was also shown by Duan 01 that nondissipative drag is responsible for an 
emergence of an interlayer Hall voltage in bilayer electron systems in the fractional quantum Hall regime. 

Basing on the previous studies one can consider nondissitive drag as a fundamental property of systems with 
macroscopic quantum coherence. For the system of uncharged bosons this effect is especially important, since the 
"drag force" plays the role similar to the role of the vector potential of magnetic field in superconductors. It opens 
new possibilities to observe the effects caused by phase coherence in such systems. One of the goals of this paper is 
to point out on this analogy. In particular, we show that the nondissipative drag effect allows to realize an entangled 
flux state in a superfluid ring with a weak Josephson link. 

Now a great attention is given to the study of ultracold alkali-metal vapours confined in magnetic and optical 
traps, where Bose-Einstein condensation of atoms was observed [TT| . Advances in technology allow to manipulate 
parameters of such systems and make ultracold atomic gases an unique object for the study of various quantum 
mechanical phenomena. 

In this paper we study the effect of nondissipative drag in a system of two quasi two-dimensional atomic Bose gases 
confined in two parallel traps. To describe such a situation we take into account that densities of atoms in the drive 
and the drag layers can be non-equal. In previous studies 0, Q only the case of two layers with equal densities of 
the particles was considered. Another important factor is the temperature. In atomic gases it is of order or higher 
than the energy of intralayer interactions. Previously, the dependence of nondissipative drag on the temperature was 
treated only qualitatively |6|, |8( • Here we study the temperature dependence quantitatively. We also evaluate the 
value of the drag for concrete mechanisms of interlayer interaction in atomic Bose gases. 



II. MODEL AND APPROACH 



The geometry of the Bose cloud can be modified significantly under a variation of a configuration of external fields 
forming a trap. When the confining potential is strongly anisotropic and the temperature as well as the chemical 
potential are smaller than the separation between the energy levels of spatial quantization in one direction the Bose 
gas can be treated as a two-dimensional one. Recently, low-dimensional atomic gases were realized experimentally 

Bose clouds of a ring shape can be created by using toroidal traps. A configuration of two toroidal traps situated 
one above another is convenient for the study of the drag effect. As it follows from the further consideration if one 
excites a circulating superflow in one trap it inevitably leads to a redistribution of this superflow between two traps 
and superfluid currents appear in both rings. 

Main features of the drag effect can be understood from the study of a system of two uniform two-dimensional Bose 
gases situated in parallel layers. The Hamiltonian of the system can be presented in the form 



/- i .2 ^ ir - i _> 



where 



is the kinetic energy, 



Ei 



I d2 ^[V^(r)]W ; (r) (2) 



Eft* = J d 2 rdV*+(r)#+(r')VHr - r')*r(r')$,(r) (3) 

is the energy of the intralayer (I — I') and interlayer (/ ^ /') interaction, and Ni — J cPr^f (r)*I';(r). Here ^ is the 
Bose field operator, I, the layer index, r, the two-dimensional radius-vector lying in the layer, and p,i, the chemical 
potentials. To be more specific we consider the case of point interaction between the atoms: Vii(r) = V22(r) = jS(r), 
Vi2(r) = V2i(r) = 7'5(r) with 7 > and |7'| < 7. Assuming the barrier between two traps is quite high we neglect 
the tunneling between the layers. 

For further analysis it is convenient to use the density and phase operator approach (see, for instance, [Tsl IT3|). 
The approach is based on the following representation for the Bose field operators 



#i(r) = exp [ifi(r) + i<pi(r)] yj pi + pi(r), 



(4) 
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®t ( r ) = Vpi exp H<M r ) - i<ft( r )] , 



(5) 



where pi and ipi are the density and phase fluctuation operators, pi = (r)^;(r)) is the c-number term of the 
density operator (one can see that it is just the density of atoms in the layer ipi(r) is the c-number term of the 
phase operator (in the approach considered the inclusion of this term in the phase operator allows to describe the 
states with nonzero average superflows). 

Substituting Eqs. into Hamiltonian JQ) and expanding it in series in powers of p\ and V^; we arrive to the 

expression 



H — Hq + Hi + H2 



In © the term 



7 2 



2- p*( v< M r )) +-^Pi~mpi 



■ 7 P1P2 



does not contain the operator part, the term 

' r s2 , 



+ J Pi + 7 P3-1 - Pi 



Pi(r) + — Pl [Vipi{r)]V<pi(r) 
m 



is linear in the phase and density fluctuation operators, and the term 

2 



H 2 = d*r(j: 



2m 



Vp ; (r) 



r + 



Kp/(v^(r)) 2 + [Vp,(r)](ft(r)Vp ( (r) + [V^(r)]pi( 
|([Vp,(r)]V0i(r) - [V^,(r)]V^(r)) + ^ [(pi(r)) 2 + (p 2 (r)) 2 ] +y^(r)^(r) 



is quadratic in V0j and /)/ operators. 

If the chemical potentials are fixed the Hamiltonian Hq is minimized under conditions 



(6) 



(7) 



(8) 



(9) 



2m 



(v^;(r)) + 7p; + 7'p3-i - Pi = 0, 



(10) 



'(piV W (r 



(11) 



Fulfillment of Eqs. (JTDJ , (fill) means that the linear in density and phase fluctuation operator term Hi in the 
Hamiltonian vanishes. One should note that, as it follows from Eq. IjlUI) . the densities of the components do not 
depend on the coordinates only when the phase gradients remain space independent as well. 

The quadratic part of the Hamiltonian determines the spectrum of elementary excitations. Hereafter we will neglect 
the higher order terms in the Hamiltonian ijfJJl . These terms describe the scattering of the quasiparticles and they can 
be omitted if the temperature is much smaller than the critical temperature (T c ~ h 2 p/m). 

III. DRAG CURRENT 

The operator of the density of the current 

f (V*+)#i - tf+V#jl , (12) 



- ih 
Jl 



2m 

rewritten in terms of the phase and density operators, has the form 



3i = — Vpi + Tl\?{vi + <Pi)W Pi + Pi- (13) 
m 
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Expanding (|13[) in powers of the density and phase fluctuation operators and neglecting the terms of order higher 
than quadratic one we obtain the following expression for the mean value of the density of the current 



^ = ~k PlVlfl + 2^ ^ V ^^ + (P* V ^» ' 



(14) 



To derive Eq. (|14|) we take into account that (<pi) = (pi) = 0. 

To compute the averages in l|14[l we rewrite the quadratic part of the Hamiltonian in terms of the operators of 
creation and annihilation of the elementary excitations. In the absence of the interlayer interaction (7' = 0) it can be 
done by the substitution 



P/( r ) 



■E 



„ikr / e k 



L[ 6j (k) + 6+(-k)] 



(15) 



(16) 



where the operators bf , bi satisfy the Bose commutation relations. Here S is the area of the system, ek = h 2 k 2 /2m 
is the spectrum of free atoms, and En~ = 'ek(ck + ^7 Pi) is the spectrum of elementary excitations at 7' = and 
Vipi = 0. 

In the case considered the substitution l|15fl . (|16fl reduces the Hamiltonian to the form 



h 2 = j: 



Ik 



fi(k)6+(k)6,(k) + -e fc ) 



J2 9k [6l"(k)6a(k) + 6i(k)6a(-k) + /i.e.] , 



(17) 



where 



£i(k) - + — kV^z, 
m 



(18) 



fffc = 7 e/c 



P1P2 



(19) 



The Hamiltonian i|17fl contains non-diagonal in Bose creation and annihilation operator terms and it can be diago- 
nalized using the u-v transformation 



&z(k) = ui a (k)a(k) + u I/9 (k)/3(k) + WiQ (k)a+(-k) + ^(k)/3+(-k) 
(see [^3]) that reduces the Hamiltonian (|17fl to the form 



(20) 



f„ (k) ( «+(k)a(k) + £)+ £p{k) (V(k)/J(k) + \ \ - , , 



(21) 



It is convenient to present the u — v coefficients and the energies of the elementary excitations as series in powers 
of gk- The u — v coefficients read as 



u la (k) u 1/3 (k) 

U2a(k) 1*2/9 (k) 



9k 

£i(k)-£ 2 (k) 



9k 



£i(k)-£ 2 (k) 
#k 



(22) 



fia(k) «i^(k)\ 

W2a(k) «2^(k) / 



where 



ffg[g 2 (k)+£ 2 (-k)] 

[£i (k)+£i (-k)] [fi(-k)-£ 2 (-k)] [£i(-k)+£ a (k)] 

~£i(-k)+£ 2 (k) 



A k =i-S 



9k 



£i(k)+£ 2 (-k) 



[£2 (k) +£ 2 (-k)] [£1 (-k) -£ 2 (-k)] [£1 (k)+£ 2 (-k)] 



2 V[^i(k)-£ 2 (k)] 2 [fi(k)+f 2 (-k)]V 



(23) 



(24) 
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B k = 1 - 



9l 



1 



1 



2 V[£i(k)-£ 2 (k)P [g 1 (_k) + f 2 (k)]V 
The spectra of the elementary excitations are found to be 

1 1 



f (k)=f 1 (k)+^ 



fi(k)-f 2 (k) fi(k)+f 2 (-k) 



(25) 



(26) 



^(k)=^(k)- fl 2 



fi(k)-f 2 (k) f x (-k)+f 2 (k) 



(27) 



One can see that the small parameter of the expansion is gk/\£i(k) — £ 2 (k)| -C 1. The last inequality takes place 
for all k, if 7' max(pi, p 2 ) <C j\pi — p 2 \- Since in most cases of interest the interlayer interaction is much smaller than 
the intralayer one, for the bilayer systems with different densities in the adjacent layers one can neglect the 0(g%) 
and higher order terms in Eqs. (EU, i|57|) . 

Using representation i|15fl . we obtain from (|14ll the following expression for the density of the current 



j, = A piV ^ + A^ k (b+(k)6 z (k)). 



(28) 



Substituting Eq. I|20|) with coefficients Q22[l. (|23(l into Eq. (|28|) . computing the averages and expanding the result in 
powers of the phase gradients we obtain the following expression for the currents: 



Jl = — [(/"si - Pdr)V(,5i + PdrVv? 2 ] , 

m 



(29) 



j 2 = — [(A2 - Pdr)^f2 + PdrV(^i] . 

m 



(30) 



Eqs. (|29|l . (|30|) are given in linear in V<p; approximation. Here the higher order in the phase gradients terms can 



be neglected if the phase gradients V<£>; are much smaller than the inverse healing lengths £ ; 



corresponds to the velocities of the superflow much smaller than the critical ones vi c ~ yfipi/'m) 
the quantities p s i and p^ r with an accuracy up to the g\ are determined by the expressions 



~ -Jmrfpi I H (that 
In Eqs. 



AU = Pi + ^ ^ 



cW, fc 



(-1)' 



cw lfe cw 2fc 



0£u 



1 



C?-E 2 /c / \ (-^lfc 



£ 2fe ) 2 
1 



(E\k — E 2 k) 2 
(-1)' 



E\k + E 2 k Eik — E 2 k 



(31) 



pdi 



1 + JV lfc + 7V 2fc Ni k - N 2k 



[Eik + E 2k f [Eik - 3 



2k I 



1 



dNik , cW 2fe 



dEu 



83 



2k 



1 



(E lk -E2 k ) 2 (E lk +E 2k y 



(32) 



Here JVjfc = [exp(Eik/T) — is the Bose distribution function. One can see that in the absence of the interlayer 
interaction (gk = 0) the value of pdr is equal to zero, and Eq. (|31l) for p s i is reduced to the well known expression for 
the density of the superfluid component at finite temperatures. If the interlayer interaction is switched on, the value 
of pdr becomes nonzero. Then, even in the absence of the phase gradient in the drag layer the superfluid current in 
this layer emerges as a response on the phase gradient in the drive layer. 

Eqs. H31|) . EH were derived under assumption of pi 7^ p 2 (and, consequently Eik 7^ E 2 k). The case pi ~ p 2 
required more rigorous consideration since in this case the mixing of the modes is strong even for a weak interlayer 
interaction. One can find that the expressions 1311) . (|32|l remain finite at pi — > p 2 : 



lim p s i 

Pl—>P2 



lim p s2 



1 ^ dNk 
S^ £k dE k 

k 



!.<? ^ 



25 



3k 



d 2 N k 
dEl 



Ek 



d 3 N k 
dEl 



(33) 
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^ ** = 45 £ I' + ^ " ^ + 3^WJ ' (34) 

where Ufc is the energy of the elementary excitations at p\ = p2 and 7' = 0. Using the exact expressions for the spectra 
and the u — v coefficients we obtain that for the case of two layers with equal densities and in the weak interlayer 
interaction limit 7' <C 7 the quantities p s and pdr are determined just Eqs. I(33() . I(34() . It allows us to conclude that 
Eqs. I|31|) . 1(32(1 are valid for an arbitrary ratio between the densities. 

Let us fist consider the case of zero temperature. We define the drag current as the current in the drag layer (e.g. 
layer 1) in the absence of the phase gradient in this layer. At T = the drag current is equal to 

idr = C dr ( — ] (7^2) — PiV(p 2 , (35) 

where 

f°° x 1/2 
C dr = / dx -3 . (36) 
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\/x + ly/x + P\j P2 {\fx + 1 + \Jx + Pl/p-^j 



The factor Cdr is an increasing function of the ratio P2I P\ (at p-ij P\ —> the factor Cdr approaches to zero, at p 2 = P\ 
it is equal to 1/12, and it approaches to 1/4 at P2I P\ — * 00). Thus, the drag current increases under increasing the 
density of the particles in the drive layer. 

At finite temperatures the drag current decreases. At small T one can use the long-wave approximation for the 
spectra Ei^)k m the temperature dependent part of Eq. I|32|l and evaluate this part analytically. It yields the following 
relation: 



jdr{T) = jdr(0) 



1 16C(3) pi / T 

Cdr P2 \ 2 7Pl 



(37) 



But, actually, this approximation is valid only at very low temperatures. The results of numerical evaluation of Eq. 
(132(1 are shown in Fig. ^ This figure demonstrates that at T ~ 27P1 a temperature decrease of the drag current is 
much slower. Basing on the results presented in Fig. ^we also conclude that the temperature reduction of the drag 
current becomes smaller under increasing the density of the particles in the drive layer. 

In a Bose cloud confined in a trap the density is nonuniform. It results in a modification of the spectrum of 
elementary excitations. We may argue that this modification reveals itself in only minor changes of the value of the 
drag. One can find that the main contribution to the sum in Eq. 132(1 comes from the excitations with the wave 
vectors of order or higher than the inverse healing lengths In systems with the healing lengths much smaller than 
the linear size of the Bose clouds the spectrum at q > S,^ 1 , ^ 1 is well described by the quasi uniform approximation. 
Therefore, in such systems the local drag current is given by the same equations l(29|) - ((32|) . as in the uniform case 
with the only modification that the densities pi, p 2 in these formulas should be understood as the functions of r. 
In particular, we predict, that for two Bose gases confined in harmonic traps having the same Thomas-Fermi radius 
(and this radius is much larger than the average healing length) a spatial distribution of the superflow in the drag 
trap T = will replicate (with a drag factor) the spatial distribution of the superflow in the drive trap. At finite 
temperatures one can expect a reduction of the drag factor near the edge of the Bose cloud, where the density is 
small. 

One can ask to which extent two-dimensionality of the system studied may influence on the results obtained. It is 
known that in 2D systems fluctuations of the phase of the order parameter are large and at nonzero temperature the 
off-diagonal one-particle density matrix (*I' + (r) ( I'(0)) goes to zero in the limit |r| — > 00. It means the absence of the 
long range order in the systems at T ^ 0. But since the asymptotic behavior of the density matrix are described by 
a power-law dependence on r (not an exponential one), at temperatures lower than the critical one (the Kosterlits- 
Touless transition temperature) the superfluid density becomes nonzero. The drag of the superflow between two 2D 
Bose-gases, considered in the paper, is connected with a finite value of the superfluid density and that is why it 
decreases under increasing the temperature. The density and phase operator approach, used in this paper, is not 
based on the existence of the Bose-Einstein condensate. Moreover, the power-law asymptotic behavior of the density 
matrix can be easily derived in this approach under accounting the thermal excitations described by the Hamiltonian 
H 2 - But at the same level of approximation we do not find any crucial influence of the two-dimensionality on the 
drag phenomena. 
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12 3 
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FIG. 1: The dependence of the drag current on the temperature at pi/ pi = 5,1,0.2 (solid, dashed and dotted curves, 
correspondingly) normalized to its value at T = 0. The dependence of jdr{0) on the ratio Pi/pi normalized to its value at 
pi = pi is shown in the inset. 



IV. THE VALUE OF THE DRAG IN ATOMIC BOSE GASES 



Let us present some estimates for the value of the drag in atomic Bose gases. For simplicity we specify the case of 
P\= Pi= p and T = 0. It is convenient to introduce the drag factor /d r = Pdr/{p — pdi) that gives the ratio between 
the currents in the drag and in the drive traps in the absence of the phase gradient in the drag trap. Taking into 
account that pdr <C p, we use the approximate expression /d r = Pdr/p for further analysis. 

The value of fd r depends on the interaction parameters 7 and 7'. The parameter 7 can be expressed through the 
dimensionless effective "scattering length" a 



7 = a. (38) 

m 

In a quasi two-dimensional trap the e ffective scattering length is connected with the 3-D scattering length a and the 
oscillator length in z direction l z — ^/hjmui z by the relation a = a/l Zl which is valid for a <C l z UM- Let us introduce 
the interlayer dimensionless effective "scattering length" a' that is connected with the interlayer interaction parameter 
7' by the relation 

i = — a'. (39) 

m 

Substituting Eqs. I|38|) . I|39|) into Eq. 134|) we obtain the drag factor in the form 
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Eq. I|40|) is valid for \a'\ <C a, but one can expect that it is approximately correct at \a'\ « a (we emphasize that in 
any case the stability condition requires the inequality \a'\ < a be satisfied). 

To evaluate the value of \a'\ we should specify a mechanism of the interlayer interaction. Let us first consider the 
interaction that corresponds to the "tail" of the Van der Waals potential: 

v.rv) - (4i) 

Here Cq is the Van der Waals constant and d is the interlayer distance. The Fourier-component of the potential (|41|) 
is 



V?r(k) = J <PrV£™(r)e«' = -^K 2 (kd), (42) 

where K2{x) is the modified Besscl function of the second kind. Taking into account that the Van der Waals interaction 
is a short-ranged one we can evaluate 7' as 7' = V^> dw (fc — > 0) = — ttCq / (2d 4 ) . It yields |&y<™r| ~ \/ tt /2C§m/ (Ah 2 d 4 ). 

This result can be obtained in a more rigorous way. Our approach is easily generalized for the case of an arbitrary 
central force interlayer interaction potential. To do this we should redefine the quantity as 



ff fc =V ia (fcW-=££-. (43) 

and substitute this definition (instead of Eq. I|19|l ) into the formulas for p s i and pdr obtained in the previous section. 
Using Eq. 1)34(1. one can present the drag factor (for T = and pi = p-i) in the following form 

J_m*_ R [°° x 2 [V 12 ( qo x)} 2 
/dr ~ 16tt 2 fi 4 5 V 2 io (x 2 + 1)5/2 ^U 

with q — y^Enpa. Substituting Eq. l|4"2")l into Eq. |H)l we find 

/dr 4(^) |y|-fVdw(d8b). (45) 

Q r 4 /•<» ,,6 

F Vd w(x) = — jf dy (1+ ^ )B/2 ir a a (sy) (46) 

describes the dependence of the drag factor on the density p. Comparing Eqs. (|46|l and l|40|) we obtain the following 
expression for the modules of the effective interlayer scattering length 

7T Ce,m 



Here the function 



Kdwl = \l 2^V^«W- ( 4? ) 



The dependence of the factor \/Fydw on the parameter dqo is shown in Fig. [21 One can see from this figure that at 
dqo <C 1 (that corresponds to the low density limit) the factor y/Fvdw m Eq. (|47|l is close to unity and we arrive to 
the expression for la^dyl gi ven above. 

Due to a short-range nature of the Van der Waals interaction the interlayer effective scattering length decreases 
quickly under increasing of d. Therefore, the interlayer distance d should be rather small to achieve an observable 
value of the drag. Using a typical value of C§ (Cq « 3 • 10~ 57 erg-cm 6 ) and taking d w lOnm , m — 87 a.u. (Rb) we 
evaluate |aydw 

I re 1CT 1 . 

The quantities l z and a can be controlled in experiments. The first one is controlled by changing the profile of the 
confining potential in z direction, and the latter one - by tuning the magnetic field to the value close to the Feshbach 
resonance field 0-01- Near this resonance the scattering length change its sign and a situation with rather small 
3D scattering length a (much smaller than the value of l z which, in its turn, has to be smaller than d/2) can be 
realized. Using this possibility one can tune the quantity a close to the value of |ay dw | and obtain the drag factor 
/,/, -- 7- ID 

Out of the resonance the typical values of 3D scattering length lie in the interval 3 -7- 5nm and for l z < d/2 and 
d pa lOnm the estimation a = a/l z is not applicable. In the ultra 2D limit (l z /a <C 1) the interaction parameter can 
be evaluated by using the formula [2(| 

Airti 2 1 



m |ln(pa 2 )| 
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FIG. 2: The density dependent factors in the effective interlayer scattering length versus the parameter dqo = \f d>npd 2 /a. 
Solid curve - \fFvdW (Van der Waals interaction), dashed curve - y/Fd-d (dipole-dipole interaction). 



For typical densities p = 10 8 -j- 10 10 cm~ 2 it yields a — 0.2 0.4 and the drag factor fd r ~ 2 3 • 10~ 3 . 

At d > lOOnm the drag caused by the Van der Waals interaction becomes negligible small. But in the last case 
the dipole-dipole interaction may give an essential contribution to the drag. Let us consider the situation where the 
dipole momenta of the atoms are aligned in a direction perpendicular to the layers. Then the interaction potential 
has the form 

^\r) = Di ^-^ /2 , (48) 

where D is the dipole momentum. The Fourier-component of the potential 148(1 reads as 

V^ d (k) = ~2TrD 2 ke- kd . (49) 

Substituting Eq. (gHJ into Eq. (@3J we obtain 

fc-s(w)'s^«(*»- (50) 



where 



poo 4 

F d - d (*) = ^ 2 l dy {l+ V y2)5/2 e-^ (51) 



One can see that Eq. (|50|l is reduced to Eq. (|40|l under definition 



l d-d\ 



D2m [* — 



(dq ). (52) 



The dependence \J 'Fd-d(dq ) is also shown in Fig. |21 In difference with the previous case the value of a' d _ d approaches 
to zero in the low density limit. But at dqo > 0.1 that corresponds to p > 10~ 2 <i~ 2 /(8wa) one can neglect the 
dependence of /<j r on the density and put the factor \jFd-d ~ 0.2. For the estimates given below we assume the 
condition dqo > 0.1 is fulfilled. 

For the magnetic dipole-dipole interaction D is the magnetic dipole momentum of the atoms. The magnetic dipoles 
can be aligned in the same direction if a constant magnetic field is applied to the system. Taking d = lOOnm, D = ps 
(the Bohr magneton) and m=87 a.u. we obtain |a^_ d | w 3 • 10 . In the case, tuning a to the |a^_J value, one can 
achieve the drag factor fd r ~ 2 • 10~ 5 

For the Bose atoms with large magnetic dipole momenta this value can be much larger. A good candidate atom is 
Cr (D = 6ps)- The possibility to realize Cr Bose-Einstein condensate is discussed in |21|. For m—52 a.u., D = 6p>B 
and d = lOOnm we evaluate |a^_J ~ 6 • 10 -3 and, consequently, the maximum drag factor f dr ~ 4 ■ 10~ 4 
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V. THE "DRAG FORCE" AS AN ANALOG OF THE VECTOR POTENTIAL 



In section IIIII we compute the drag current directly. The same results can be obtained from the analysis of the 
dependence of the free energy of the system on the phase gradients. The free energy of the system can be found from 
the common thermodynamic relation 



F = H + E zero + T Y, ln 

\=a,(3 k 



1 — exp 



T J 

Here the quantity Hq given by Eq. Q is the classical energy of the system, and 



(53) 



E 



= \ E £^(k)-e fe ] (54) 



A=q,/3 k 

is the energy of the zero-point fluctuations. 

Substituting the spectra (|26(l , (|27(l into Eq. H53J) and expanding the final expression in powers of the phase gradients 
we find the following expression for the free energy 

f h 2 

F = F Q + / d 2 r— [p s i(V(pi) 2 + p s2 {Vip 2 ) 2 - Pdr(V>i - Vt^) 2 ] + higher order terms, (55) 

where Fo does not depend on the phase gradients and the quantities p B i, pdr are determined by the expressions 
(|31[) . (|32|) . One can see that the answer (|29[) - l|32l) obtained in Sec. IIIII by another method can also be found from Eq. 
(|55|l using the relation 

1 9F /KM 

hs o(\/ipi) 

The relation (|55l) is more instructive in a sense that it demonstrates an analogy between the drag effect in superfluids 
and the exciting of a supercurrent by an external magnetic field in superconductors. To illustrate this analogy let us 
consider two ring-shape traps and fix the phase gradient in the drive trap (trap 2 in further notations). Then the free 
energy as the function of the phase gradient in the drag trap (trap 1) can be presented in the form 

F = const + — /5 s i($ + $dr) , (57) 

mn 

where R is the radius of the ring, w, its width, p s i = p s \ — pdr, 

$ = -!-</ eflVVi (58) 



2tt j c 

(here C is a contour around the ring) is the winding number for the phase tpi and 

$dr = / rflV^ 2 , (59) 

Psl 27T J c 

is the winding number for the phase tp 2 times the drag factor. In deriving (|57|l we, for simplicity, neglect the dependence 
of densities on the coordinate inside the traps. 

Since the value of $ should be integer the minimum of the free energy at |$dr| < 1/2 is reached for $ = 0. In 
this case the phase gradient in the drag trap is equal to zero and the superfluid current in the drag trap flows in the 
same direction as in the drive trap. If |$dr| > 1/2 the free energy reaches its minimum at nonzero $ and the phase 
gradient is induced in the drag trap. Then together with the drag current the counterflow current appears in the drag 
trap (depending on the value of <E>dr the total current in this trap can be parallel as well as antiparallel to the current 
in the drag trap). Just the same situation takes place in a superconducting ring with nonzero flux of magnetic field 
inside the ring. Thus, in two-ring Bose systems the quantity <I>dr plays the same role as a flux of an external magnetic 
field (measured in units of flux quanta) in superconducting circuits. 

To realize this situation experimentally one should create a circulating supcrflow in the drive trap. It can be done 
by elliptic rotating deformation of this trap. The rotation can be switched of when a superflow be created. The value 
of drag current can be found from measurement of the angular momentum of the drag trap. At present a number 
methods for measuring this quantity has been realized experimentally |22 | -|25 | . The methods are based on the study 
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of dynamics of collective excitations, on the investigation of interference phemonena under hyperfine state transitions 
and on the observation of the dynamics of expansion of the Bose cloud. 

To extend the analogy with superconductors let us consider the case where the drag trap of a ring geometry contains 
a weak Josephson link. Then the free energy as the function of the phase shift Atp on the link reads as 



where Ej is the Josephson energy. At Ej > (Ej) c = h 2 wp sl / (2nRm) and |$dr| — 1/^) 3/2, . . . the dependence F(A<p) 
has two degenerate minima. If Ej / {Ej) c — 1^1 these minima are very shallow and one can expect that two quantum 
states with different phase shifts (and with the superfluid currents flowing in opposite directions) will be entangled. 
It is the same regime that is required for implementing the superconducting Josephson flux (persistent current) qubit 
[2fil |. While in alkali- metal Bose gases the drag factor is rather small and even in the most favorable conditions the 
maximum value can be reached is of order 10~ 2 -^10~ 3 (see Sec. II V|l . the case |$dr| ~ 1/2 can be realized in ring-shape 
traps of large radiuses (10 2 ~ 10 3 fim). 

In conclusion, we would like to mention another systems in which the effects described in this paper may take place. 
It is excitonic or electron-hole Bose liquids in electron bilayers. In these systems electron-hole pairs with components 
belonging to adjacent layers may form a superfluid state. For the first time the effect was predicted in [27| 28], and 
recently it was confirmed experimentally |29j. The superfluid drag effect may emerge in two parallel bilayers (the 
four-layer system). In the four-layer system the intralayer (in the same bilayer) and interlayer (between bilayers) 
interactions are of the same order: both of them arc determined by the dipole-dipole mechanism. In such a case the 
dipolc momentum of the pair is large. Therefore, one can expect that nondissipative drag in these systems will be 
rather strong. 

This work is supported by the INTAS grant No 01-2344. 
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